AD-A035  196  GEORGIA  INST  OF  TECH  ATLANTA  SCHOOL  OF  ELECTRICAL  EN-- ETC  F/G  12/2 

ON  AN  OPERATOR  THEORY  OF  LINEAR  SYSTEMS  WITH  PURE  AND  DISTRIBUT— ETC (U) 
1975  E W KAMEN  DA-AR0-D-31-124-73-G171 

UNCLASSIFIED  AR0-U554.6-M  NL 


3-77 


SECURITY  Cl  ASSI  E|C  ATION  OF  THIS  PAGE  (Wh*t,  Dtlt  Eritorod; 


Jr'  READ  INSTRUCTIONS 

UEFOKE  COMPLETING  yphM 

3.  RECIPIENT'S  CATALOG  NUMBER 


REPORT  DOCUMENTATION  PAGE 


12  GOVT  ACCESSION  NO. 


5.  TYPE  OF  REPORT  d PERIOD  COVERED 


TiTLE  (mid  Subtitle) 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AuTHORi'iy 


CONTRACT  OR  GRANT  NUMBERS; 


PERFORMING  organization  NAME  AND  AODRESS 


10.  PROGRAM  ELEMENT,  PROJECT.  TASK 


AREA  A WORK  UNIT  NUMBERS 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


13-  numberofpage: 


4 MONITORING  NGEnCY^AME  a kDORESSOtdlTterent  from  Controlling  Office) ! 15.  SECURITY  CLASS,  (ot  thie  report) 


Unclassified 


I5«.  DECLASSI  FI  CATION  'DOWN  GRADING 
SCHEDULE 


It.  O.STRiHUTtON  ST  AT  EMENT  (o!  this  H»portj 


Approved  for  public  release;  distribution  unlimited 


17.  OiSTRifciuTiON  STATEMENT  (of  the  mbetrect  entered  In  Block  20,  If  different  from 


d Supplementary  notes 

The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position,  unless  so  designated  by  other  authorized 

documents. 


19.  KEY  WORDS  (Continue  on  reveree  etde  If  neceeeery  end  identify  by  block  number) 


Representation  theory 
Delay  equations 


Convolution 


20.  ABSTRACT  (Continue  on  reveree  etde  If  neceeeery  end  Identity  by  block  number) 


A representation  theory  based  on  convolution 
operations  is  developed  for  a large  class  ot* 
linear  systems  containing  pun  and  distributed 
delays  in  state  and  control.  in  terms  of  this 
framework  a necessarv  and  sufficient  condition 
and  a sufficient  condition  are  given  for  func- 
tional (null)  controllabil ity.  The  conditions 
involve  the  generation  of  modules  defined  over 
a convolution  ring  of  functions. 


Unclassified  ' u 3 ^ w 

SECURITY  Cl  ASSlElCATlrtN  r»F  TmiS  PAftF  fWhen  nmte  Knteted) 


WA5-2 


in  Proceedings  of  the  1975  IEEE  Conference  on  Decision  and 
Control,  Houston,  Texas,  pp.  77-80,  1975. 

ON  AN  OPERATOR  THEORY  OF  LINEAR  SYSTEMS  WITH  PURE  AND  DISTRIBUTED  DELAYS* 


E . W . Kamen 

School  of  Electrical  Engineering 
Georgia  Institute  of  Technology 
Atlanta,  Georgia  30332 


A representation  theory  based  on  convolution 
operations  is  developed  for  a large  class  of 
linear  systems  containing  pure  and  distributed 
delays  in  state  and  control.  In  terms  of  this 
framework  a necessary  and  sufficient  condition 
and  a sufficient  condition  are  given  for  func- 
tional (null)  controllability.  The  conditions 
involve  the  generation  of  modules  defined  over 
a convolution  ring  of  functions. 


1.  Introduction 

In  many  control  problems  the  systems  under 
consideration  contain  pure  and  distributed  time 
delays  in  state  and  control  (examples  are  given 
by  MANITIUS  [1]).  Such  systems  are  usually  re- 
ferred to  as  hereditary  systems  since  the  rate  of 
change  of  the  present  state  depends  on  past  values 
of  the  state  and  control  or  input. 

In  this  paper  we  consider  the  class  of  linear 
systems  given  by  a first-order  functional  differen- 
tial equation  of  the  form 

r°  r 

x(t)  = I A(9)x(t  + 9)de  + F0x(t)  + T FjXU  - 

-c  i=l 


+ I B(9)u(t+9)d9  + G u(t)+  lG  u(t-b  ) 
J-d  i-1  1 1 


where  c,d  and  the  a^,b^  are  positive  real  numbers, 
the  F^  (G.)  are  nxn  (nym)  matrices  over  the 

reals  R,  A(9)  (resp.  B(9))  is  a nxn  (nxm)  matrix 
of  (Lebesgue)  measurable  and  integrable  functions 
on  [-c,0]  ([-d,0]),  x(t)  6 Rn  is  the  "instantaneous 
state,”  and  u(t)  € R"1  is  the  input. 


DELFOUR-McCALLA-MITTER  [5]  (see  also  the  survey 
by  MANITIUS  [6]). 

In  contrast  to  existing  methods,  our  approach 
to  the  study  of  (1)  is  based  on  an  algebraic 
setting  defined  in  terms  of  convolution  operators. 
More  precisely,  in  the  next  section  it  is  shown 
that  (1)  can  be  written  in  the  form 

x(t)  - (F*x)(t)  + (G  * u)  (t)  (2) 

where  * denotes  convolution  and  F and  G are 
matrices  whose  elements  belong  to  a convolution 
ring  of  functions  and  impulses  (Dirac  distribu- 
tions). The  convolution  representation  (2)  is  a 
special  case  of  the  time-domain  operator  framework 
developed  by  KAMEN  *7]. 

In  the  latter  part  of  the  paper  the 
representation  (2)  is  applied  to  the  problem  of 
driving  initial  functions  to  the  zero  function  in 
finite  time  (functional  null  controllability). 

New  algebraic  criteria  for  controllability  are 
given  in  terms  of  modules  defined  over  the  convo- 
lution rings. 

2.  Representation  by  Convolution  Operators 
Xoc 

Let  L denote  the  space  of  all  real-valued 
Lebesgue  measurable  functions  f(t)  that  are 

locally  integrable,  i.e.  I |f(t)|dt  < <*  for  any 

K loc 

compact  subset  K of  R.  Let  L denote  the  sub- 
1 oc 

space  of  L consisting  of  all  functions  with 
support  bounded  on  the  left.  It  is  easily 

verified  that  L*°  is  a ring  with  pointwise 
addition  and  with  convolution  defined  by 


m 

(g*f)(t)  - [ g(9) f (t  - 9)d9 


Systems  of  the  form  (1)  have  been  studied 
using  mainly  functional-analytical  methods  applied 
to  a state  space  setting  defined  in  terms  of  the 
product  space  Rn x L^ (-h,0; Rn)  where  h “ maxfc.a^}. 

In  particular,  numerous  results  on  controllability 
and  optimal  feedback  control  can  be  found  in  the 
work  of  DELFOUR-MITTER  [2,3],  DELFOUR  [4],  and 
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Note  that  if  the  support  of  g is  contained  in 
[0,a',  a > 0,  then 

„a 

(g*f)(t)  = g(9)f(t-6)d9 

J o 

and  by  a change  of  variables,  we  have 


3.  Initial  Conditions 

In  this  section  we  show  how  to  "handle" 
nonzero  initial  conditions  in  solving  the  operator 
equation  x « F*x  + G*u.  This  result  will  be 
utilized  in  the  next  section  to  study  functional 
controllability. 


(g*  f)(t)  = f g(0)f(t+  0)d9  (3) 

-a 

where  g(9)  = g(-9). 

Letting  5^  denote  the  unit  impulse  (Dirac 
distribution)  concentrated  at  the  point  [a),  define 

q . 

J = it+  Z b 6 :f£L°C,  a ,b  € R, 

L i-1  1 ai  1 1 

q “ positive  integer^-  . (4) 

The  addition  in  the  definition  of  J can  be  taken 
in  a formal  sense  or  it  can  be  viewed  as  addition 
in  some  space  of  distributions.  The  set  J is  an 
loc 

overring  of  L+  with  the  convolution  operation 
given  by 

'vf  + Ibi6a.)*(g  + rci4d.)  = 

f * g + Z (c^f (t  - d^  + bjgCt  - ai)) 

+ ^(^bjVa.+d.)  • <5> 

j i J J i 

Note  that  5 is  the  identity  element  of  the  ring  J. 

Finally,  let  L£  denote  the  space  of  all  real- 
valued measurable  functions  f(t)  defined  a.  e.  on 
R with  support  bounded  on  the  left,  such  that 

I | f (t ) | dc  < ®.  The  space  L is  a subring  of  L_7°c. 

■I*  + + 

Via  the  above  constructions  we  can  now 
characterize  (1)  in  terms  of  convolution  operations. 

Theorem  1 : With  u€  (L£)m  and  x€  Jn,  the  system 

equation  (1)  can  be  written  in  the  form 
x<t)  = (F*x)(t)  + (G*u)(t)  where  F (G)  is  a 
n*n  (nxm)  matrix  over  Jc  = (at  J: supp  a is 

compact  and  contained  in  [0,®)]. 


Suppose  that  we  are  given 
x(t)  « (F  * x) (t)  + (G  * u) (t) 


(6) 


with  F and  G defined  over  J as  in  Theorem  1. 

c 

Now  let  a = max  {r  £ supp  F)  and 
b « max  (t  6 supp  G}.  Then  to  solve  (6)  for 
t > 0 we  need  to  specify  x(t)  on  the  interval 
[-a,0]  and  u(t)  on  the  interval  [-b,0].  We  assume 
the  following  initial  conditions 


x(0) 


X e r 

o 


Oil 


x(t)  - 0(t),  t€[-a,0),  5)  € ^L^_a 
U(t)  - uo(t),  t([-b,o],  uoe  (LJ.b>0^) 

where  ll  d-  « {f€L£:supp  fc[c,d]}. 


(7) 


Given  f 6 L+,  define 


£l[c,d](t) 


f(t),  t € [c,d] 


0,  otherwise 


By  the  results  of  DELFOUR  [4],  with  initial  data 
/I  \ra 

(7)  and  u|,  . 6IL.  . ; , the  convolution 

' (o,«)  t°.">  / 

equation  (6)  has  a unique  solution  x(t)  with 


!(o,«)  €(L[o%))  • 


(l}  . 'j  with  initial 

\ [o,®>) 


Theorem  2:  Given  u| , , € 

I (o,«) 

data  (7),  the  solution  of  (6)  for  t > 0 is  equal 
to  the  solution  for  t > 0 of 

x(t)  = (F  * x) (t)  + (G  * u) (t)  + v(t) 

with  Initial  data  equal  to  zero,  where 


v = x 5 +(F*0  + G*u)l.  . € 

o o o (o,®)  c 


Proof : Applying  (3)  and  (5)  to  (1),  we  have  that 

* * £ Ft8  >x+  (£  + f Ci6  )*u 

i*o  l i=o  i 

where  A(t)  = A(-t)  and  8(t)  = B(-t). 

By  definition  of  A,B  and  the  F^.G^  the  coefficient 

matrices  of  x and  u are  over  J . ■ 

c 


Proof:  Clearly  xq6o  sets  up  the  initial  value 

x € Rn.  Since  x “ os  + x|  . . and  u = u|  , ,+u 

o 1 (o,«)  1 (o,®)  o 

with  x ■ 0,  we  have  that 
o 

F*  x + G*  u = F*  x| , . + G * u|  . +F*^  + G*u 

(o,®)  1 (°,®)  o 


Hence 

(F*x  + G*u)!(o>.)  = F*xI(o>-)+G*u|(o>.) 

+ <F^  + C*Uo)!(o,«)  . • 
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As  a consequence  of  Theorem  2,  operational 
methods,  such  as  those  developed  in  [7],  can  be 
used  to  solve  operational  differential  equations 
of  the  form  (6)  with  nonzero  initial  conditions. 

We  shall  now  apply  this  result  to  the  study  of 

controllability. 

4.  Controllability 

Let  L1  denote  the  set  of  all  f € L*  such  that 
c + 

supp  f is  compact  and  contained  in  [o,«).  With  the 
induced  operations  L*  is  a convolution  ring  con- 
tained in  the  ring  J . 

Definition:  The  system  given  by  (6)  is  said  to 

be  (null)  controllable  if  for  any  initial  con- 

1 ra 

ditlon  (7),  there  exists  a control  u € (L  ) such 

c 

that  the  solution  of  (6)  is  zero  for  all  t > h, 
some  h > 0. 

In  terms  of  the  following  constructions  we  derive 
a necessary  and  sufficient  condition  for  control- 
lability. 

Let  pn  denote  the  nt*1  derivative  of  6 in  the 

n o 

sense  of  distributions.  The  element  p belongs  to 
J2T',  the  convolution  ring  of  Schwartz  distributions 

on  R with  support  bounded  on  the  left  (see  [7]). 

Given  f € L+  , the  derivative  of  f in  the  sense 

of  distributions  is  equal  to  p * f (with  the  con- 
volution p*f  carried  out  in.£p.  Therefore  (6) 

can  be  expressed  completely  in  terms  of  convolu- 
tion operations: 


denote  the  L^-submodule  of  J [p]n  generated  over 
1 c c 
Lc  by  Yx.  v2 . That  is,  <Vj V >L1 

i c 

is  the  set  of  all  sums  T-V.y, , tt,  €L  . 

ll  i c 

Theorem  3:  The  system  given  by  (6)  is  controll- 
able if  and  only  if  Jn  is  contained  in  the  L*- 

c 1 c 
module  generated  over  Lc  by  the  columns  of  the 

nxn  matrix  (pi  - F)  and  the  nXra  matrix  C;  that 
is, 


J"  C <(pl  - F),G)l1 


(9) 


l.m 


Proof:  Suppose  that  (9)  holds  and  the  initial 

condition  (7)  is  given.  Then  there  exist  ar€  (L*) 
and  S€  (L^)n  such  that 

(pi  - F)*8  - G*ar  - xq60+  (F*0  + G*uq)|(o  m)  6 j". 

Therefore  by  (8)  and  Theorem  2,  g is  the  solution 
of  (6)  with  the  given  initial  data  and  with  input 
u)  m j “ or.  Hence  the  system  is  controllable. 

Conversely,  suppose  that  the  system  is 
controllable. 

Define  e » (0  0 ...  1 0 ...  0)TR€  Rn,  i-  l,2,...n. 

4 ,th  , 
i-i  place 

By  Theorem  2,  e^6Q  sets  up  the  initial  value 

x = e, . Then  since  the  system  is  controllable, 

® ^ 1 m 

by  (8)  and  Theorem  2 there  exist  u,€(L  ) and 
, i n l c 

x,t(L  ) such  that 
i c 


(pi  - F)  *x  =*  G*u  (8) 

where  I is  the  nxn  identity  matrix. 

Now  let  J [p]  denote  the  set  of  all  finite 
sums  E»j*P  where  oo  €.1^.  With  the  standard 
operations,  Jc[p]  is  a convolution  ring  of  dis- 
tributions (contained  in  XF') . The  rings  J and 
1 + c 

L (cJ  ) are  subrings  of  J [pi. 
c c c J 

Letting  Jc[p]n  denote  the  space  of  n-element 
column  vectors  over  Jc[p],  we  have  that  Jc[p]n  Is 
a free  finite  module  over  J^fp]  with  componentwise 
addition  and  with  multiplication  defined  by 

"(gj,  e2,  ....  en)TR' («*H1(t»*s2 "*en>TR 

where  rr,  Pj€Jc[p],  and  TR  denotes  the  transpose. 

Since  L^  is  a subrtng  of  .1  [pj,  by  restricting 

c C J 

the  multiplication  operation  to  L we  have  that 

1 

I .Pj  ip  al=>  a (nonfiuilei  L^-module  Given 
V*,»  Vj ^ let  <Y]»  V2* 


(pi  - F)  *xi  - G*u1  - eii0>  1“  1,2 n. 

But  e.fi  , e.6  , ...,  e 6 is  a basis  of  jn  as  a 
1 o 2 o no  c 

J^-module,  so  that  given  yij^,  there  exist 

v,  £ J such  that  v » £ y.  (e. S ).  Hence 
i c 'i  i o 

y - (Pi  - F)*(Svi*xi)  - C*  (Eyi*u1).  (10) 

Then  since  o*f€L1  for  any  o€J  , fCL1,  it 
c 1 c c’ 

follows  from  (10)  that  the  columns  of  (pi  - F)  and 

G generate  Jn  over  L.  ■ 
c c 

Corollary;  The  system  (6)  is  controllable  if  and 
only  if 

ei&o€  <(pl  - F),  C>  1#  1*1,  2 n 

LC 

where  e *(00...  1 0 ...  0) € Rn 


[_tth 


place 


In  the  finite-dimensional  case  where  F and  G 
are  over  R6q,  it  is  well  known  (KALMAN  [8])  that 

•-he  system  given  by  (6)  is  controllable  if  and 
only  if  the  rnnk  of  the  nynm  matrix 

(C,FG F"'^)  is  equal  to  n.  If  F ard  G are 
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viewed  as  matrices  over  R (rather  than  Rfi  ) , this 

o 

is  equivalent  to  requiring  that 

<G,FG,...,Fn‘lG>R  - Rn  . (11) 

Let  N be  a fixed  subring  of  J with  6 €N. 

c o 

In  view  of  (11)  it  Is  reasonable  to  ask  If  the 
condition 

(C,FC Fn'lC>„  » Nn  (12) 

is  necessary  and  sufficient  for  controllability 
when  F and  G are  over  N.  The  answer  is  that  (12) 
is  sufficient  but  not  necessary.  The  proof  of 
sufficiency  will  be  given  in  an  expanded  version 
of  this  paper.  The  simple  example  below  shows 
that  (12)  is  not  necessary. 

Example : Suppose  that  x(t)  = x(t)  + u(t  - a),  a>0. 

In  this  case  F **  6 and  G * 6 are  over  the  subring 
o a 

RC s ] » fZb,4,  :b  U]  c J . 
k aJ  1 i la  1 ‘ c 

Since  the  inverse  of  J is  4 ,(4  + N for 

a -a’  a N 

S =■  R[6  ] or  N ■ J . 

'•a-  c 

However  the  system  is  controllable  since  there 

exist  Q,  8 € L*  such  that  (p  - 4 ) * 0+ 6 *8“  4 . 

c o a o 

For  Instance  we  could  take 
eC  , 0 t t * a 

o(t)  = /_ef_(t  - b)  ,a*  t S b 
J b-a 

0,  otherwise 


“ — ( t + a-b-1),  0*t*b-a 
b-a  ’ 


0,  otherwise 


where  b is  any  fixed  real  number  greater  than  a. 

The  function -8  given  above  is  a control  for  the 
initial  condition  x(0)  - 1. 

5.  Further  Applications 

In  addition  to  function  controllability, 
there  exist  many  problem  areas  that  can  be  investi- 
gated using  the  convolution  representat ion  (6). 

For  example,  the  results  of  SONTAG  [9]  can  be 
carried  over  to  this  framework  giving  a necessary 
and  sufficient  algebraic  criterion  for  Euclidean 
reachability. 

A particularly  interesting  topic  is  the  study 
of  feedback  control  systems  with  u • K*x, 

K » mXn  matrix  over  a subring  N of  J^.  Practical 

problems  Include  the  development  of  algebraic 
procedures  for  the  design  of  K's  to  achieve 
stabilization  or  pole  allocation.  Results  on  this 
are  already  available  in  the  case  that  F and  G are 


over  a subring  N that  is  a principal  ideal  domain 
(see  MORSE  [10]  and  SONTAG  [9]). 
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